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q \ Abstract 

£sQ ■ The statistical expectation values of the temperature fluctuations of cosmic microwave 

background (CMB) are assumed to be preserved under rotations of the sky. We inves- 
tigate the statistical isotropy of the CMB anisotropy maps recently measured by the 
Wilkinson Microwave Anisotropy Probe (WMAP) using bipolar spherical harmonic 
power spectrum proposed in Hajian & Souradeep 2003. The Bipolar Power Spectrum 

(N ' 

(BiPS) is estimated for the full sky CMB anisotropy maps of the first year WMAP 
data. The method allows us to isolate regions in multipole space and study each re- 
gion independently. This search shows no evidence for violation of statistical isotropy 
in the first-year WMAP data on angular scales larger than that corresponding to 

■ 1 Introduction 

(N 

■ In standard cosmology, CMB anisotropy signal is expected to be statistically isotropic, i.e., statistical 
expectation values of the temperature fluctuations AT(q) are preserved under rotations of the sky. In par- 
ticular, the angular correlation function C(q, q') = (AT(q)AT(q')) is rotationally invariant for Gaussian 
fields. In spherical harmonic space, where AT(q) = J2i m a imYim((}) the condition of statistical isotropy 
(SI) translates to a diagonal (ai m a*, m ,) = Ci6u>5 mm > where C;, the widely used angular power spectrum 
of CMB anisotropy. SI CMB sky is essential for C; to be a complete description of (Gaussian) CMB 
anisotropy and hence an adequate measure for comparing with models. Hence, it is crucial to be able 
to determine from the observed CMB sky whether it is a realization of a statistically isotropic process, 
or not. The detection of statistical isotropy (SI) violations in the CMB signal can have exciting and 
far-reaching implication for cosmology. For example, a generic consequence of cosmic topology is the 

^ ■ breaking of statistical isotropy in characteristic patterns determined by the photon geodesic structure of 

the manifold as probed by the CMB photons traveling to us from the surface of last scattering over a 
distance comparable to the cosmic horizon, Rh- On the other hand, SI violation could also arise from 
foreground contamination, non-cosmological signals and be artifacts of observational technique. 

The first-year Wilkinson Microwave Anisotropy Probe ( WMAP) observations are consistent with pre- 
dictions of the concordance ACDM model with scale- invariant and adiabatic fluctuations which have been 
generated during the inflationary epoch [Hinshaw et al. 20031 |Kogut et al. 2003| |Spergel et al. 2003| 
|Page et al lol^ |FeIrTs et al., 2003] . After the first year of WMAP data, the SI of the CMB anisotropy 
(i.e. rotational invariance of n-point correlations) has attracted considerable attention. Tantalizing 
evidence of SI breakdown (albeit, in very different guises) has mounted in the WMAP first year sky 
maps, using a variety of different statistics. It was pointed out that the suppression of power in the 
quadrupole and octopole are aligned |Tegmark et al. 2004| . Further "multipole- vector" directions associ- 
ated with these multipoles (and some other low multipoles as well) appear to be anomalously correlated 
|Copi et al. 20041 ISchwarz et al. 2004| . There are indications of asymmetry in the power spectrum at 
low multipoles in opposite hemispheres [Eriksen et al. 2004al IHansen et al. 20041 |Naselsky et al. 2004 . 
Possibly related, are the results of tests of Gaussianity that show asymmetry in the amplitude of the mea- 
sured genus amplitude (at about 2 to 3a significance) between the north and south galactic hemispheres 
[Park 2004IIEriksen et al. 2004bllEriksen et al. 2004cj . Analysis of the distribution of extrema in WMAP 
sky maps has indicated non-gaussianity, and to some extent, violation of SI [Larson fc WandelT 2004 . 
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However, what is missing is a common, well denned, mathematical language to quantify SI (as distinct 
from non Gaussianity) and the ability to ascribe statistical significance to the anomalies unambiguously. 

Since the observed CMB sky is a single realization of the underlying correlation, the detection of SI 
violation or correlation patterns pose a great observational challenge. For statistically isotropic CMB 
sky, the correlation function 

C(nx,n 2 ) = C(ni -n 2 ) = ^2 / dKC{Kh u TZn 2 ), (1) 

where IZh denotes the direction obtained under the action of a rotation 1Z on n, and dTZ is a volume 
element of the three-dimensional rotation group. The invariance of the underlying st^isticjijmder rotation 
allows the estimation of C(hi ■ 77-2) using the average of the temperature product AT(n)AT(n') between 
all pairs of pixels with the angular separation 9. In the absence of statistical isotropy, C(h, n') is estimated 
by a single product AT(n)AT(n') and hence is poorly determined from a single realization. Although 
it is not possible to estimate each element of the full correlation function C(n,n'), some measures of 
sjjatistical^ anisotropy of the CMB map can be estimated through suitably weighted angular averages of 
AT(h)AT(h'). The angular averaging procedure should be such that the measure involves averaging over 
sufficient number of independent 'measurements', but should ensure that the averaging does not erase 
all the signature of statistical anisotropy (as would happen in eq. or eq. J2U). Recently, we proposed 
the Bipolar Power spectrum (BiPS) m (£ = 1, 2, 3, . . .) of the CMB map as a statistical tool of detecting 
and measuring departure from SI |Hajian fe Souradeep 2003b| |Souradeep fc Hajian 2003| and reviewed 
in this article in sec. The non-zero value of the BiPS spectrum imply the break down of statistical 
isotropy 

STATISTICAL ISOTROPY n e = W ^ 0. (2) 

The BiPS is sensitive to structures and patterns in the underlying total two-point correlation function 
|Hajian fe Souradeep 2003rj| |Souradeep fc Hajian 2003| . The BiPS is particularly sensitive to real space 
correlation patterns (preferred directions, etc.) on characteristic angular scales. In harmonic space, the 
BiPS at multipole £ sums power in off-diagonal elements of the covariance matrix, (ai rn ai> m i}, in the 
same way that the 'angular momentum' addition of states Im, I'm' have non-zero overlap with a state 
with angular momentum \l — l'\ < £ < I + I'. Signatures, like ai m and a;+ nm being correlated over a 
significant range I are ideal targets for BiPS. These are typical of SI violation due to cosmic topology 
and the predicted BiPS in these models have a strong spectral signature in the bipolar multipole £ 
space |Hajian fc Souradeep 2003a] . The orientation independence of BiPS is an advantage since one can 
obtain constraints on cosmic topology that do not depend on the unknown specific orientation of the 
pattern (e.g., preferred directions). 

The results of WMAP are a milestone in CMB anisotropy measurements since it combines high angular 
resolution, high sensitivity, with 'full' sky coverage allowed by a space mission. The frequency coverage 
allows for WMAP CMB sky maps to be foreground cleaned up to I ~ 100 |Tegmark et al. 2004| . The 
CMB anisotropy map based on the WMAP data are ideal for testing for statistical isotropy. 



2 Sources of Statistical Isotropy violation 

An observed map of CMB anisotropy, AT° bs , contains the true CMB temperature fluctuations, AT,, 
convolved with the beam and buried into noise and foreground contaminations. The observed map AT 
is related to the true map through this relation 

AT obs = Y^B qj AT 3 + N t , (3) 

3 

in which B is a matrix that contains the information about the beam smoothing effect and n is the 
contribution from instrumental noise and foreground contamination. Hence, the observed map is a 
realization of a Gaussian process with covariance C = C T + C N + C res where C T is the theoretical 
covariance of the CMB temperature fluctuations, C N is the noise covariance matrix and C res is the 
covariance of residuals of foregrounds. Breakdown of statistical isotropy C(n,n') ^ C(n ■ n') can occur 
in any of these parts of the correlation function. Broadly, these effects may be divided into two kinds: 
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• Theoretical signals: These effects are theoretically motivated and are intrinsic to the true CMB 
sky, AT. We discuss two examples of these effects, i.e. non-trivial cosmic topology and primordial 
magnetic fields, in the next subsections. 

• Observational artifacts: In an ideally cleaned CMB map, the true CMB temperature fluctuations 
are completely extracted from the observed map. But this is not always true. Sometimes there are 
some artifacts (related to B or N) left in the cleaned map which may in principle violate the SI. 
These effects are explained in section l2~3l 

2.1 Cosmic Topology & Ultra-large scale structure 

The cosmic microwave background anisotropy is currently the most promising observational probe of 
the global spatial structure of the universe on length scales near to and even somewhat beyond the 
'horizon' scale (~ cHq 1 ). Figure ^ depicts a prevalent view within the concept of inflation, that this 
relatively smooth Hubble volume that we observe is perhaps a tiny patch of an extremely inhomogeneous 
and complex spatial manifold. The complexity could involve non-trivial topology (multiple connectivity) 
on these ultra-large scales. Within a general program to address the observability of such a diverse 
global structure, a more well defined and tractable path would be to restrict oneself to spaces of uniform 
curvature (locally homogeneous and isotropic FRW models) but with non-trivial topology; in particular, 
compact spaces which have additional theoretical. 

The question of size and the shape of our universe are very old problems studied earlier [Ellis 197T1 
ISokolov fc Shvartsman 19741 ^ |Lachieze-Rey fc Luminet~9 95 . With remarkable improvements in cos- 
mological observations, in particular the CMB anisotropy measurements, these questions have received 
considerable attention over the past few years |de Oliveira-Costa et al. 19* 96 , Starkman 1998 , Le vin et al. 19981 
|Bond,Pogosyan fc Souradeep 1998,2000|lde Qlivei~ Costa et al. "200llDineen et al. 2004l|Copi et al. 2004| . 
Although a multiply connected universe sounds non-trivial, but there are theoretical motivations |Linde" 2004 
ILevin 2002j to favor a spatially compact universe. One possibility to have a compact flat universe is the 
consideration of multiply connected (topologically nontrivial) spaces. The oldest way of searching for 
global structure of the universe is by identifying ghost images of local galaxies and clusters or quasars 
at higher redshifts |Lachieze-Rey fc Luminet 1995| . This method can probe the topology of the universe 
only on scales substantially smaller than the apparent radius of the observable universe. Another method 
to search for the shape of the universe is through the effect on the cosmic density perturbation fields. 
For compact topologies, the two main effects on the CMB are: (1) the breaking of statistical isotropy in 
characteristic patterns determined by the photon geodesic structure of the manifold and (2) an infrared 
cutoff in the power spectrum of perturbations imposed by the finite spatial extent. More generally, in 
a universe with non-trivial global spatial topology, the multiple connectivity of the space could lead to 
observable characteristic angular correlation patterns in the CMB anisotropy. 

Over the past few years, many independent methods have been proposed to search for evidence of a 
finite universe in CMB maps. These methods can be classified in three main groups. 

• Using the angular power spectrum of CMB anisotropies to probe the topology of the Universe. 
The angular power spectrum, however, is inadequate to characterize the peculiar form of the 
anisotropy manifest in small universes of this type. Since nontrivial topology breaks down SI, there 
is more information in a map of temperature fluctuations than just the angular power spectrum 
[Levm et al. 19981 |Bond,Pogosyan fc Souradeep 1998,2000| |Hajian fc Souradeep 2003a| . 

• The second class of methods are direct methods that rely on multiple imaging (or strong correlation 
features) of the CMB sky. The most well known methods among these methods are S-map statistics 
jde Oliveira-Costa et al. 19961 Ide Oliveira-Costa et al. 207)3] and the search for circles-in-the-sky 
[C"orni sh,Spcr gel fc Starkman 1998| . 

• Third class of methods are indirect probes which deal with the correlation patterns of the CMB 
anisotropy field by using an appropriate combination of coefficients of the harmonic expansion of 
the field [Dineen et al. 20041 |Donoghue et al. 2004| |Hajian fc Souradeep 2003a]|Copi et al. 2004| . 
The Bipolar power spectrum (BiPS) method is one of the strategies in this class. 
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Figure 1: A cartoon depicting a prevalent view within the inflationary paradigm. The observable universe 
corresponds to a small patch of a very complicated manifold that has been blown to cosmological scales 
during an inflationary epoch. Ultra-large scale structure could be observable if the the size of this patch 
is not much smaller that the scales of inhomogeneity and non-trivial topology. 

The correlation patterns in CMB that lead to violation of SI implies that imply (a/ m a^ m ) has off- 
diagonal elements. Figure El taken from |Bond,Pogosyan fc Souradeep 1 998,2000| shows the off-diagonal 
elements in the CMB correlation for two compact universe models. BiPS gathers together the power in 
the off-diagonal elements of (ai m a* m } as shown in Fig. 

Using the fact that statistical isotropy is violated in compact spaces one could use the bipolar power 
spectrum as a probe to detect the topology of the universe. A simple example of is the BiPS signature 
of a non-trivial topology can be given for a T 3 universe, where the correlation function is given by 

C(q, q') = L- 3 p *( fc n) e-"( £ ' n -9-V n -9'), (4) 

in which, n is 3-tuple of integers (in order to avoid confusion, we use q to represent the direction instead 
of n), the small parameter < 1 is the physical distance to the SLS along q in units of L/2 (more 
generally, L/2 where L = (L1L2L3) 1 / 3 ) and L is the size of the Dirichlet domain (DD). When e is a small 
constant, the leading order terms in the correlation function eq. (@J can be readily obtained in power 
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Figure 2: The figure taken from |Bond,P ogosya n~fe Souradeep 1998,2000| illustrates the non-diagonal 
nature of the expectation values of ag m pair products when the CMB anisotropy violates SI in two 
model compact universe. The radical violation in the model on the left corresponds to a small compact 
universe where CMB photons have traversed across multiple times. The model on the left with mild 
violation of SI corresponds to a universe of size comparable to the observable horizon. For more details, 
see |Bond,Pogosyan fc Souradeep 1998,20001 



series expansion in powers of e. For the lowest wave numbers |n| 2 = 1 in a cuboid torus 

C(q,q') « 2^P <& (2 7 r/ J L l )cos(7re/3 t A< Zl ) 



(5) 



Co 



l-e 2 |A# + 3e 4 £(A^ 



where Aqi are the components of Aq = q — q' along the three axes of the torus and Pi = L/ Li. From 
this, the non-zero Ki can be analytically computed to be 
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Ki has the information of the relative size of the Dirichlct domain and one can use it to constrain 
the topology of the universe. A detailed study of the BiPS signature of cosmic topology is given in 
|Hajian fc Souradeep 20 03a . These prediction allow us to constrain cosmic topology using the BiPS 
measured in the observed CMB maps |Hajian et al. 2004| . 



2.2 Primordial Magnetic Fields 

Cosmological magnetic field, generated during an early epoch of inflation Ratra 1992, Bamba et al. 2004 , 
can generate CMB anisotropies [Durrer et al. 1998| . The presence of a preferred direction due to a 
homogeneous magnetic field background leads to non-zero off-diagonal elements in the covariance ma- 
trix [Chen et al. 2004] . This induces correlations between a; +l rra and a/_i jTO multipole coefficients of the 
CMB temperature anisotropy field in the following manner 

(ai m a*, m ,) = 5 m , m >[6tj>Ci + (6t + ij'-i + Si-ij'+iDt)], (7) 
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where Di is the power spectrum of off-diagonal elements of the covariance matrix. For a Harrison-Peebles- 
Yu-Zel'dovich scale-invariant spectrum, D\ behaves as l~ 2 . More precisely, it is given by 

D l = 4 x l(T 16 r 2 (-^) 4 . (8) 

This clearly violates the statistical isotropy and gives rise to a non-zero BiPS predictions for magnetic 
fields. This open the way to use BiPS analysis on CMB maps to constrain or measure primordial cosmic 
magnetic fields |Hajian et al. 2004b| . 



2.3 Observational Artifacts 

Foregrounds and observational artifacts (such as non-circular beam, incomplete/non-uniform sky coverage 
and anisotropic noise) would also manifest themselves as violations of SI. 

• Anisotropic noise : The CMB temperature measured by an instrument is a linear sum of the cosmo- 
logical signal as well as instrumental noise. The two point correlation function then has two parts, 
one part comes from the signal and the other one comes from the noise 



C(n 1; h 2 ) = C s ( ni , h 2 ) + C N '(h lt h 2 ). (9) 

Both signal and noise should be statistically isotropic to have a statistically isotropic CMB map. 
So even for a statistically isotropic signal, if the noise fails to be statistically isotropic the resultant 
map will turn out to be anisotropic. The noise matrix can fail to be statistically isotropic due to 
non-uniform coverage. Also if the noise is correlated between different pixels the noise matrix could 
be statistically anisotropic. A simple example of this is the diagonal (but anisotropic) noise given 
by the following correlation 

C N (h,n') = <y 2 (n)5 nfl ,. (10) 
This noise clearly violates the SI and will lead to a non-zero BiPS given by 

e 

*t= Yl \f*™\ 2 > ( n ) 

m=—t 

where fi m are spherical harmonic transform of the noise, fi m = J dflfiY^ (n)a (n). 

• The effect of non-circular beam : In practice when we deal with data, it is necessary to take into 
account the instrumental response. The instrumental response is nothing but the beams width and 
the form of the beam and can be taken into account by defining a beam profile function B(h, h'). 
Here n denotes the direction to the center of the beam and n' denotes the direction of the incoming 
photon. The temperature measured by the instrument is given by 



AT (h) = / AT(h')B(h, h')dn w . (12) 



Using this relation to calculate the correlation function C(h\, h-z) = (AT (hi) AT (112)) one would 
get 



C(h x ,h 2 ) = I dVL lV J <m A »(AT(h')AT(n"))B(hi, h')B(h 2 , h") (13) 
(Kin 1 I dQfi''C(n',n'')B(fii, n')B(n 2 , h"). 



Only for a circular beam where B(n, h') = B(n ■ n'), the correlation function is statistically 
isotropic, C(hi, h 2 ) = C(h\ ■ h 2 ). Breakdown of SI is obvious since even Ci get mixed for a 
non-circular beam, Ci = Y^v Ml'Cv |Mitra et al. 2004] . Non-circularity of the beam in CMB 
anisotropy experiments is becoming increasingly important as experiments go for higher resolution 
measurements at higher sensitivity. 
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Mask effects : Many experiments map only a part of the sky. Even in the best case, contamination 
by galactic foreground residuals make parts of the sky unusable. The incomplete sky or mask effect 
is another source of breakdown of SI. But, this effect can be readily modeled out. The effect of a 
general mask on the temperature field is as follows 

AT masked (h) = AT(n)W{h), (14) 

where W(fi) is the mask function. One can cut different parts of the sky by choosing appropriate 
mask functions. Masked a; m coefficients can be computed from the masked temperature field, 



aZ sked = J AT masfced (n)Y^(n)dfi ft (15) 
= X! a 'i m i J Y hm 1 {n)Y* m (n)W{h)dD, fl . 



l\ TTll 



Where ai l7ni are spherical harmonic transforms of the original temperature field. We can expand 
W(n) in spherical harmonics as well 

W{n) = £)u;j m Y, m (n), ( 16 ) 
lm 

and after substituting this into eq. I)15|l it is seen that the masked ai m is given by the effect of a 
kernel K 1 ^ 1 on original a tm [Prunet et al. 2004| 

masked \ ^ rylitni /i t\ 

a lm = 2^ a *i™i^m ■ ( 17 ) 
Iimi 

The kernel contains the information of our mask function and is defined by 

K\T = ^^r^jY hmi (h)Y hr ^(h)Y^(h)dSh (18) 

(2f! + l)(2b + l) l0 lm 

4ir(2l + 1) O (i0/20 Ly iimii 2 m 2 ■ 

The covariance matrix of a masked sky will no longer have the diagonal form because of the action 
of the kernel 



1-211)2 




I masked masked 
\ a lm a l'm' 



*} = (a Wmi ^)K\™I$$ (19) 
= C^i'Mm'Ktr^C' 1 

li ,mi 

This clearly violates the SI and results a non-zero BiPS for masked CMB skies. In the next section 
we apply a galactic mask to ILC map and show that signature of this mask on BiPS is a rising tail 
at low ^ (£< 20). 

Residuals from foreground removal : Besides the cosmological signal and instrumental noise, a CMB 
map also contains foreground emission such as galactic emission, etc. The foreground is usually 
modeled out using spectral information. However, residuals from foreground subtractions in the 
CMB map will violate SI. Interestingly, BiPS does sense the difference between maps with grossly 
different emphasis on the galactic foreground. As shown in |Hajian fc Souradeep 2005b| the BiPS 
of a Wiener filtered map shows a signal very similar to that of a galactic cut sky. This can be 
understood if one writes the effect of the Wiener filter as a weight on the 'contaminated' galactic 
regions of the map. 

AT w (fi) = AT(h)(l + W(h)). (20) 

This explains the similarity between a cut sky and a Wiener filtered map. The effect of foregrounds 
on BiPS still needs to be studied more carefully. 
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3 The Bipolar Power Spectrum (BiPS) 



Two point correlation of the CMB anisotropy is given by ensemble average, but there is only one observ- 
able CMB sky. Hence, the ensemble average is meaningless unless the CMB sky is SI, when the two point 
correlation function C(9) can be well estimated as in eq. ()21|l by the average product of temperature 
fluctuations over ahjpairs of directions h\ and fi2 whose angular separation is 9. In particular, for CMB 
temperature map AT(nj) defined on a discrete set of points on celestial sphere (pixels) hi (i — 1, . . . , N p ) 

N P 

C(6) = ^2 AT(ni)AT(n 3 -)<J(cos0-nf -hj), (21) 

is an estimator of the correlation function C(9) of an underlying SI statistics. If the statistical isotropy is 
violated the estimate of the correlation function from a sky map given by a single temperature product 

C{hi,h 2 ) = AT(n!)AT(n 2 ) (22) 

is poorly determined. 

Although it is not possible to estimate each element of the full correlation function C(ni,fi2), some 
measures of statistical isotropy of the CMB map can be estimated through suitably weighted angular 
averages of AT(fii)AT(fi2)- The angular averaging procedure should be such that the measure involves 
averaging over sufficient number of independent measurements , but should ensure that the averaging 
does not erase all the signature of statistical anisotropy. Another important desirable property is that 
measure be independent of the overall orientation of the sky. Based on these considerations, we have 
proposed a set of measures of statistical isotropy |Hajian fc Souradeep 20035| 

K e = (21 + l) 2 J dQ ni J dfl n2 [JL J dTZ X £ (n)C(TZh u Kh 2 )} 2 . (23) 

In the above expression, C{TZh\, IZh^) is the two point correlation at 1Zh\ and lZfi2 which are the 
coordinates of the two pixels h\ and fi2 after rotating the coordinate system through an angle u> where 
(0 < u> < 7r) about the axis n(0, $). The direction of this rotation axis n is defined by the polar angles 
where (0 < O < 7r) and <&, where (0 < $ < 2ir). \ l 1S the trace of the finite rotation matrix in the 
^M-representation 

X \K) = J] flU(R), (24) 

M—-£ 

which is called the characteristic function, or the character of the irreducible representation of rank I. 
It is invariant under rotations of the coordinate systems. Explicit forms of y^(TV) are simple when 1Z is 
specified by ui, 0, $, then x^(^) 1S completely determined by the rotation angle ui and it is independent 
of the rotation axis n(0, $), 

X e (TZ) = x\") (25) 
sin[(2£ + l)w/2] 
sin[w/2] ' 

And finally dlZ in ea. (|23|l is the volume element of the three-dimensional rotation group and is given by 

dU = 4 sin 2 - dio sin 9 d® tf$ . (26) 

For a statistically isotropic model C(ni, fi2) is invariant under rotation, and therefore C(TZh\, lZfi2) = 
C(ni, fi2) and the orthonormality of x^(w), we will recover the condition for SI, 

if = k°*«). (27) 

Real-space representation of BiPS is very suitable for analytical computation of BiPS for theoretical 
models where we know the analytical expression for the two point correlation of the model, such as 
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theoretical models in |Hajian fc Souradeep 20 03a . On the other hand, the harmonic representation of 
BiPS that we describe next allows computationally rapid methods for BiPS estimation from a given CMB 
map. 

Two point correlation of CMB anisotropics, C(ni, 112), is a two point function on S 2 x S 2 , and hence 
can be expanded as 

C(n 1 ,h 2 )= A uA Y h(ni)®Y l2 (n 2 )} m , (28) 

h,l 2 ,L,M 

where A\¥ are coefficients of the expansion (here after BipoSH coefficients) and {Yi 1 (fii) ® Yi 2 (n 2 )}^M 
are the bipolar spherical harmonics which transform as a spherical harmonic with I, M with respect to 
rotations Varshal ovich et al. 1988) given by 

{5 / !l (ni)®4(«2)}»/= C hLhn l2 Y hmi (n2)Y hm2 (h 2 ), (29) 

mim2 

in which Cf 1 ^^ m are Clebsch-Gordan coefficients. We can inverse-transform C(fii, h 2 ) to get the Aj^f 
by multiplying both sides of ea. l|25|l by {^(ni) <g) Yi' 2 (n 2 )}g, M , an d integrating over all angles, then the 
orthonormality of bipolar harmonics implies that 

4Y 2 = / dfifti / d(l n2 C(ni, n 2 ) {Y h (n x ) ® Y h (n 2 )}| M . (30) 



The above expression and the fact that C(h\, n 2 ) is symmetric under the exchange of h\ and h 2 lead to 
the following symmetries of A\H^ 



Af M = { _ l)(h+l2 - L)Kl (31) 
£M r 

II °i,2k+li 



A e n M = Ai M 5, 2k+1 , k = 0,1,2,.... 



The Bipolar Spherical Harmonic (BipoSH) coefficients, Ayf , are linear combinations of off-diagonal 
elements of the covariance matrix, 

A lll 2 ~ { a hmi a l 2 m 2 )^^) m2( ^lim 1 l2-rn2- (32) 

ra\ra 2 

This means that A\H completely represent the information of the covariance matrix. Fig. [21 shows how 
Af¥ and Af¥ combine the elements of the covariance matrix. When SI holds, the covariance matrix is 
diagonal and hence 

A e u ¥ = (-1)^(2/ + 1) 1 / 2 ^, S tQ S M o, (33) 

BipoSH expansion is the most general representation of the two point correlation functions of CMB 
anisotropy. The well known angular power spectrum, Ci is a subspace of BipoSH coefficients correspond- 
ing to the A®® that represent the statistically isotropic part of a general correlation function. When SI 
holds, A\f or equivalently C; have all the information of the field. But if the SI breaks down, Aff are 
not adequate for describing the field, and one needs to take the other terms into account. This simply 
means that the And when the statistical isotropy holds, these coefficients will reduce to the well-known 
angular power spectrum of CMB anisotropy. 

It is impossible to measure all Af¥ individually because of cosmic variance. Combining BipoSH 
coefficients into Bipolar Power Spectrum reduces the cosmic variance 3 . BiPS of CMB anisotropy is 
defined as a convenient contraction of the BipoSH coefficients 



= £ KT>o. (34) 



IJ'M 



3 

variance 



This is similar to combining a; m to construct the angular power spectrum, C\ = £m l a im| 2 i to reduce the cosmic 
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Figure 3: BipoSH coefficients are linear combinations of elements of the covariance matrix. Here Afft' 1 
(left) and Ajft 1 (right) are plotted to show how BiPS covers the off-diagonal elements of the covariance 
matrix in harmonic space. 

The BiPS, which can be shown that is equivalent to the one in ea. l|23(l . has interesting properties. It 
is orientation independent and is invariant under rotations of the sky. For models in which statistical 
isotropy is valid, BipoSH coefficients are given by eq. (|33|l . And results in a null BiPS, i.e. K£ — for 
every positive £, 

Kt = K S ea . (35) 

3.1 Unbiased Estimator of BiPS 

An estimator for measuring BipoSH coefficients from a given CMB map is 

iff - J2 VWW^ai m a llm , Cf%, m , , (36) 

mm' 

where Wi is the Legendre transform of the window function. The above estimator is a linear combination 
of Ci and hence is unbiased. An unbiased estimator of BiPS is given by 



Y.W ( 37 ) 



ll'M 



where the bias for the BiPS is defined as 93^ = (ki) — K£ is equal to 

Q5^ = ^ W h W h E [{ a i 1 m 1 ai im ' 1 )(ai 2m2 a hm ' 2 ) + (ah^h^O ^h^hrnO] 

l\,l2 mi ,m! T712 ,T7l2 

EnlM nlM fnn\ 

^l 1 m 1 l 2 m2^lim[l 2 m 2 • \°° ) 



X 

M 



The above expression for 25<> is obtained by assuming Gaussian statistics of the temperature fluc- 
tuations. The procedure is very similar to computing cosmic variance (which is discussed in the next 
section), but much simpler. However, we can not measure the ensemble average in the above expression 
and as a result, elements of the covariance matrix (obtained from a single map) are poorly determined 
due to the cosmic variance. The best we can do is to compute the bias for the SI component of a map 

i+h 

^ i = (kf) sl = (2e+l)J2 E C h C l2 W h W l2 [l + (-l) e S hl2 ] . (39) 
h h = \i-h\ 
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Note , the estimator kg is unbiased, only for SI correlation, i. e. , (kg) = 0. Consequently, for SI correlation, 
the measured kg will be consistent with zero within the error bars given by a SI |Hajian fc Souradeep 2 003b . 
We simulated 1000 SI CMB maps and computed BiPS for them using different filters. The average BiPS 
of SI maps is an estimation of the bias which can be compared to our analytical estimation. The left panel 
of Fig. 0] shows that the theoretical bias (computed from average C/) match the numerical estimations 
of average m of the 1000 realizations of the SI maps. 

It is important to note that bias cannot be correctly subtracted for non-SI maps. Non-zero kg estimated 
from a non-SI map will have contribution from the non-SI terms in full bias given in eq. (SB) . It is not 
inconceivable that for strong SI violation, 25^ over-corrects for the bias leading to negative values of kg. 
What is important is whether measured kt differs from zero at a statistically significant level. 




Figure 4: Left: Analytical bias for a Gaussian window function with W / i G (40) computed from the average 
Ci from 1000 realizations of a SI CMB map compared with ^ K j eahzatton ^ (the average ki from 1000 
realizations). This shows that the theoretical bias is a very good estimation of the bias for a statistically 
isotropic map. Right: The cosmic error, a(ng), obtained using 1000 independent realizations of CMB 
(full) sky map matches the analytical results shown by dotted curve with triangles . This shows a 
much better fit to the theoretical cosmic variance compared to what was obtained for 100 realizations 
Ha j i an fc Souradeep 2003b| 



3.2 Cosmic Variance of BiPS 

A crucial point is how well one can hope to estimate the BiPS given the single observed sky. This is 
limited by the Cosmic variance of the BiPS estimator defined as 

(T 2 = < kj > - < kg > 2 (40) 

It is possible to obtain an analytic expression variance of k using the Gaussianity of AT. Looking back 
at the ea. (|23() we can see, we will have to calculate the eighth moment of the field 

(AT(f ll )AT(n 2 )AT(n 3 )AT(fM)AT(n 5 )AT(n 6 )AT(n 7 )AT(fi 8 )). (41) 

Assuming Gaussianity of the field we can rewrite the eight point correlation in terms of two point 
correlations. One can write a simple code to do that 4 . This will give us (8 — 1)!! = 7x5x3= 105 terms. 
These 105 terms consist of terms like: 

(AT(n 1 )AT(n 2 ))(AT(n 3 )AT(n 4 ))(AT(n 5 )AT(n 6 )}(AT(rV)AT(n 8 )), (42) 

and all other permutations of them. On the other hand (kg) has a 4 point correlation in it which can 
also be expanded versus two point correlation functions. If we form (kg) — (kg) 2 , only 96 terms will be 

4 F90 software implementing this is available from the authors upon request 
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left which are in the following form 

{^^-f J dn l ---d^ i J dnj dR'x t (R)xi{R')C{hu R'n A )C{n 2 , fl / n 3 )C(i?ni,n4)C( J Rn2,n 3 ) (43) 

and all other permutations. As described in detail in our paper |Hajian fc Souradeep 2005b| , it is possible 
to simplify and group together the 96 terms and obtain a compact expression as 



l:2l>( 
i+h 



2 ^irf + { ~ 1)e{2e + 1} + (1 + 2 (- 1 )') F » 



E E AC h c h w h w h K 2(: + !) + Kh] 



16 i 



(21 + 1) 2 
2/i + l 

-tf E 



(2£+l) 



i+h 



2h + 1 



E 

=ii-iii 

E cla 2 wlw h 



(44) 



Numerical computation of cr^ is fast. But the challenge is to compute Clebsch-Gordan coefficients 
for large quantum numbers. We use drc3j subroutine of netlib 5 in order to compute the Clebsch-Gordan 
coefficients in our codes. Again we can check the accuracy of our analytical estimation of cosmic variance 
by comparing it against the standard deviation of BiPS of 1000 simulations of SI CMB sky. The result 
is shown the right panel of in Fig. 0] and shows a very good agreement between the two. 



4 Results of BiPS analysis of WMAP CMB maps 

We carry out measurement of the BiPS, on the following CMB anisotropy maps 

A) a foreground cleaned map (denoted as TOH') |Tegmark et al. 2004| , 

B) the Internal Linear Combination map (denoted as 'ILC in the figures) [Bennett et al. 2003) . and 

C) a customized linear combination of the QVW maps of WMAP with a galactic cut (denoted as 
'CSSK'). 

Also for comparison, we measure the BiPS of 

D) a Wiener filtered map of WMAP data (denoted as 'Wiener') |Tegmark et al. 2004| , and 

E) the ILC map with a 10° cut around the equator (denoted as 'Gal. cut.'). 

Angular power spectra of these maps are shown in Fig. [SJ The best fit theoretical power spectrum 
from the WMAP analysis 6 |Spergel et al. 2003| is plotted on the same figure. Ci from observed maps 
are consistent with the theoretical curve, C[ (except for the lowest multipoles). The bias and cosmic 



variance of BiPS depend on the total SI angular power spectrum of the signal and noise C\ = Cf 



■c?. 

However, we have restricted our analysis to I < 60 where the errors in the WMAP power spectrum is 
dominated by cosmic variance. It is conceivable that the SI violation is limited to particular range of 
angular scales. Hence, multipole space windows that weigh down the contribution from the SI region of 
multipole space will enhance the signal relative to cosmic error, a si . We use simple filter functions in I 
space to isolate different ranges of angular scales; a low pass, Gaussian filter 



(/ s ) = exp(-(Z + l/2) 2 /(^ + l/2) 2 ) 



(45) 

5 http: / / www.netlib.org/slatec / sre / 

6 Based on an LCDM model with a scale-dependent (running) spectral index which best fits the dataset comprised of 
WMAP, CBI and ACBAR CMB data combined with 2dF and Ly-a data 
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/ 

Figure 5: Top: Cg of the two WMAP CMB anisotropy maps. The red, magenta and green curves 
correspond to map A, B and C, respectively. The black line is a 'best fit' WMAP theoretical Ct used 
for simulating SI maps. Blue dots are the average C/ recovered from 1000 realizations. Bottom: These 
plots show the window functions used. The dashed curves with increasing I coverage are 'low-pass' 
filter, W^(l s ), with l s — 4,18,40, respectively. The solid lines are 'band-pass' filter W^kJs) with 
(l s J t ) = (13, 2), (30, 5), (30, 20), (45, 20), respectively. 

that cuts off power on small angular scales (< l/l s ) and a band pass filter, 

Wf(l t ,l s ) = [2(1 - J ((Z + l/2)/(l t + 1/2)))] exp(-(Z + l/2) 2 /(/ s + 1/2) 2 ) (46) 

that retains power within a range of multipoles set by It and l s . The windows are normalized such that 
Y,t(l + V 2 )/W + = 1> i-e-, unit rms for unit flat band power d = 1/(1(1 + 1)). The window 

functions used in our work are plotted in figure We use the Cf to generate 1000 simulations of 
the SI CMB maps. az m 's are generated up to an l ma x of 1024 (corresponding to HEALPix resolution 
Nside = 512). These are then multiplied by the window functions W^(l s ) and Wf (lt,l s )- We compute 
the BiPS for each realization. Fig[S]shows that the average power spectrum obtained from the simulation 
matches the theoretical power spectrum, Cf, used to generate the realizations. We use Cf to analytically 
compute bias and cosmic variance estimation for k^. This allows us to rapidly compute BiPS with la 
error bars for different theoretical Cj ' . 

We use the estimator given in ea. (|36|l to measure BiPS for the given CMB maps. We compute the 
BiPS for all window functions shown in Fig El Results for one these windows are plotted in Figs. In 
the low-/ regime, where we have kept the low multipoles, BiPS for all three given maps are consistent 
with zero. But in the intermediate-/ regime (Fig. [SJl, although BiPS of ILC and TOH maps are well 
consistent with zero, the CSSK map shows a rising tail in BiPS due to the galactic mask. To confirm 
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Figure 6: Measured BiPS for maps A, B and C filtered with a window with l s — 30, It = 20. This is 
to check the statistical isotropy of the WMAP in the modest 20 < I < 40 range in the multipole space 
where certain anomalies have been reported. ILC with a 10-degree-cut (top) has the same BiPS as map 
C (/ s = 30, k = 20) which explains that the raising tail of CSSK map is because of the mask. 

it, we compute the BiPS for the ILC map with a 10-degree cut around the galactic plane (filtered with 
the same window function). The result is shown on the top panel of Fig. HJ Another interesting effect 
is seen when we apply a Wf ^(20,45) filter, where Wiener filtered map has a non zero BiPS very similar 
to that of CSSK but weaker. The reason is that Wiener filter takes out more modes from regions with 
more foregrounds since these are inconsistent with the theoretical model. As a result, a Wiener filtered 
map at W*(20, 45) filter has a BiPS similar to a cut sky map. The fact that Wiener map has less power 
at the Galactic plane can even be seen by eye! Hence using different filters allows us to uncover different 
types of violation of SI in a CMB map. In our analysis we have used a set of filters which enables us to 
probe SI breakdown on angular scales I < 60. 

The BiPS measured from 1000 simulated SI realizations of Cj is used to estimate the probability 
distribution functions (PDF), p{kg). A sample of the PDF for two windows is shown in Fig. [TJ Measured 
values of BiPS for ILC, TOH and CSSK maps are plotted on the same plot. BiPS for ILC and TOH 
maps are located very close to the peak of the PDF. We compute the individual probabilities of the map 
being SI for each of the measured ke- This probability is obtained by integrating the PDF beyond the 
measured kg. To be precise, we compute 

poo 

P(ki\Cf) = P(n e >ki)= I dKtp(Kt), k e > 0, (47) 

J ke 
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Figure 7: Probability distribution function for m to K5 constructed from 1000 realizations. The left 
panel shows the PDF for the maps filtered with Wf (20,30) (left panel) and Wf(l s = 40) (right panel). 
The latter is more skewed, which explains the apparent ~ la shift in the Ki values for Wj (l s — 40) at 
low £. The green, magenta and red (circular, pentagonal and rectangular) points represent ILC, CSSK 
and TOH maps, respectively. The smooth solid curves are Gaussian approximations. 



The probabilities obtained are shown in Figs. 0] and |U for W s '(20,30), W G (40) and W, G (4). The proba- 
bilities for the Wf (20, 30) window function are greater than 0.25 and the minimum probability at ~ 0.05 
occurs at K4 for W^ G (40). The reason for systematically lower SI probabilities for Wi S (20,30) as com- 
pared to Wi (40) is simply due to lower cosmic variance of the former. The contribution to the cosmic 
variance of BiPS is dominated by the low spherical harmonic multipoles. Filters that suppress the ai m 
at low multipoles have a lower cosmic variance. 

It is important to note that the above probability is a conditional probability of measured kg being SI 
given the theoretical spectrum Cj (used to estimate the bias). A final probability emerges as the Bayesian 
chain product with the probability of the theoretical Cf used given data. Hence, small difference in these 
conditional probabilities for the two maps are perhaps not necessarily significant. Since the BiPS is close 
to zero, the computation of a probability marginalized over the Cj may be possible using Gaussian (or, 
improved) approximation to the PDF of k^. 

The important role played by the choice of the theoretical model for the BiPS measurement is shown 
for a Wi that retains power in the lowest multipoles, I = 2 and I = 3. Assuming C ; T , there are hints of non- 
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SI detections in the low £'s (top- left panel of Fig.|5J|. We also compute the BiPS using a Cj for a model 
that accounts for suppressed quadrupole and octopole in the WMAP data Shaf ieloo fc Souradeep 2004 . 
The mild detections of a non zero BiPS vanish for this case (top-right panel of Fig. El - The corresponding 
SI probabilities for the two choices of Cj are shown in the lower panels. 

5 Discussion and Conclusion 

The SI of the CMB anisotropy has been under scrutiny after the release of the first year of WMAP 
data. We use the BiPS which is sensitive to structures and patterns in the underlying total two-point 
correlation function as a statistical tool of searching for departures from SI. We carry out a BiPS analysis 
of WMAP full sky maps. We find no strong evidence for SI violation in the WMAP CMB anisotropy maps 
considered here. We have verified that our null results are consistent with measurements on simulated 
SI maps. The BiPS measurement reported here is a Bayesian estimate of the conditional probability of 
SI (for each n e of the BiPS) given an underlying theoretical spectrum Cj . We point out that the excess 
power in the Cj with respect to the measured C; from WMAP at the lowest multipoles tends to indicate 
mild deviations from SI. BiPS measurements are shown to be consistent with SI assuming an alternate 
model Cj that is consistent with suppressed power on low multipoles. Note that it is possible to band 
together ki measurements to tighten the error bars further. The full sky maps and the restriction to low 
I < 60 (where instrumental noise is sub-dominant) permits the use of our analytical bias subtraction and 
error estimates. The excellent match with the results from numerical simulations is a strong verification 
of the numerical technique. This is an important check before using Monte-Carlo simulations in future 
work for computing BiPS from CMB anisotropy sky maps with a galactic mask and non uniform noise 
matrix. 

There are strong theoretical motivations for hunting for SI violation in the CMB anisotropy. The possi- 
bility of non-trivial cosmic topology is a theoretically well motivated possibility that has also been observa- 
tionally targeted jEllis 19711 |Lachieze-Rey fc Luminet 1995|ILevin 20021 iLinde 2004] , The breakdown of 
statistical homogeneity and isotropy of cosmic perturbations is a generic feature of ultra large scale struc- 
ture of the cosmos, in particular, of non trivial cosmic topology |Bond,Pogosyan fc Souradeep 1998,2000| . 
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Figure 9: Top; Figure compares the measured values of ki for maps A and B filtered to retain 
power only on the lowest multipoles, I = 2 and I = 3 assuming the WMAP theoretical spectrum 
WMAPbf (left) and a model spectrum that matches the suppressed power at the lowest multipoles 
IShafieloo fe Souradeep 2004) . The non zero K£ 'detections' assuming the WMAP theoretical spectrum 
become consistent with zero for a Cf that has power suppressed at low multipoles. Bottom: The cor- 
responding SI probability assuming the WMAP theoretical spectrum, Cf (left) and a model spectrum 
that matches the suppressed power at the lowest multipoles (right). 



The underlying correlation patterns in the CMB anisotropy in a multiply connected universe is related 
to the symmetry of the Dirichlet domain. The BiPS expected in flat, toroidal models of the uni- 
verse has been computed and shown to be related to the principle directions in the Dirichlet domain 
|Hajian fc Souradeep 20 03a . As a tool for constraining cosmic topology, the BiPS has the advantage of 
being independent of the overall orientation of the Dirichlet domain with respect to the sky. Hence, the 
null result of BiPS can have important implication for cosmic topology. This approach complements direct 
search for signature of cosmic topology |Cornis h,Spergel & Star kman 1998|lde Oliveira-Costa et al. 1996| 
and our results are consistent with the absence of the matched circles and the null S-map test of the 
WMAP CMB maps [Cornish et al. 2003llde Oliveira-Costa et al. 2003) . Full Bayesian likelihood compar- 
ison to the data of specific cosmic topology models is another approach that has applied to COBE-DMR 
data |Bond,Pogosyan fc Souradeep 1 998,2000 ! . Work is in progress to carry out similar analysis on the 
large angle WMAP data. We defer to future publication, detailed analyzes and constraints on cosmic 
topology using null BiPS measurements, and the comparison to the results from complementary ap- 
proaches. There are also other theoretical scenarios that predict breakdown of SI that can be probed 
using BiPS, e.g., primordial cosmological magnetic fields [Durrer et al. 1998llChen et al. 2004 . 
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The null BiPS results also has implications for the observation and data analysis techniques used to 
create the CMB anisotropy maps. Observational artifacts such as non-circular beam, inhomogeneous 
noise correlation, residual stripping patterns, etc. are potential sources of SI breakdown. Our null BiPS 
results confirm that these artifacts do not significantly contribute to the maps studied here. Foreground 
residuals can also be sources of SI breakdown. The extent to which BiPS probes foreground residuals 
is yet to be fully studied and explored. We do not see any significant effect of the residual foregrounds 
in ILC and the TOH maps as it was mentioned by [Eriksen et al. 2004cj . This can not be necessarily 
called a discrepancy between the two results unless we know what should have been seen in the BiPS. 
The question is if the signal is strong enough and whether the effect smeared out in bipolar multipolc 
space within our angular /-space window. On the other hand, the very fact that BiPS does show a strong 
signal for the Wiener filtered map, mean that at some level BiPS is sensitive to galactic residuals. 

In summary, we study the Bipolar power spectrum (BiPS) of CMB which is a promising measure of 
SI. We find null measurements of the BiPS for a selection of full sky CMB anisotropy maps based on 
the first year of WMAP data. Our results rule out radical violation of statistical isotropy in the CMB 
anisotropy measured by WMAP. 
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